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We study different realisations of the first order deconfinement phase transition inside a compact 
star by comparing the Gibbs and Maxwell construction for the mixed phase. The hadronic sector is 
described within the relativistic mean field model including hyperons. The quark sector is described 
by the MIT Bag model. We find that these two realisations lead to very different star properties, 
in particular, the composition of the stellar matter. We also find that for the Maxwell construction 
there is a sharp discontinuity in the baryon density and the electron chemical potential. We argue 
that a sharp jump in the elctron chemical potential should lead to the redistribution of electrons 
and formation of strong electric fields around the discontinuity surface. 
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The heavy ion experiments at RHIC, LHC and FAIR are designed to study strongly interacting matter 



> 

' under extreme conditions of high temperature and/or high baryon density. These experiments are also 

expected to shed light on the properties of a new phase of strongly interacting matter, the Quark Gluon 
Plasma (QGP) and on the nature of the deconfinement phase transition. On the other hand, compact 
stars, due to their large central densities, serve as a natural laboratory to study the properties of the 
strongly interacting matter at high densities and small temperature, in particular, the possibility of a 
ON . deconfining phase transition. 

At finite temperature and zero baryon density, numerical studies on lattice are believed to provide 
reliable results for the physics of the deconfinement transition [l|. In this case the lattice calculations 
predict that the deconfinement happens via a smooth crossover transition Q at a temperature ^170 
- 200 MeV |3|. However, the studies at finite baryon densities on the lattice are very difficult. Some 
progress has been made in recent years in extending calculations to finite quark chemical potentials but 
they do not provide reliable results yet 0, 0] ■ There is an indication of a critical point at a rather small 
quark chemical potential fi q ta 100 MeV, with a first order transition for larger /j, q Qj. 

The most recent interest in the study of first order deconfining phase transition is related to the nature 
of the mixed phase (MP) Q. Especially, the role of the screened Colulomb potential and interface effects 
in the MP were studied by several authors (see e.g. refs. A detailed study employing the 

Wigner Seitz cell approach [7| suggests that the MP behaves more in accordance with that for Maxwell 
Construction (MC) rather than the Gibbs Construction (GC). Earlier similar questions were addressed 
in connection with the " pasta" phases associated with the liquid-gas phase transition [12j . 

In this work we present a comparative study of compact stars with the MC and GC of the mixed 
phase. For the hadronic phase we use a Relativistic Mean Field (RMF) model of the Walecka type 
Besides nucleons, this model contains hyperons as well as hyperon-hyperon interaction. For the quark 



sector we employ the MIT Bag model 1J| which has been used previously for the description of strange 



quark stars (see a recent review [15|)- We construct the MP for the deconfinement phase transition and 
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obtain the equations of state (EOSs) for the both cases. The next section is devoted to the description 
of these models and corresponding EOSs. In section 3 we use these EOSs to describe properties of the 
compact stars. In section 4 we study the jump in the electron chemical potential at the interface between 
the two phases and estimate an induced electric field. In the last section we summarise our results. 



II. PROPERTIES OF MATTER IN COMPACT STARS 
A. Hadron phase 

At low densities the relevant degrees of freedom are hadrons. To describe the hadronic phase we use a 
non-linear version of the RMF model. In this model the baryons interact with mean meson fields. The 
variant that we use here is known as the TM1 model [TgI ] . 

The Lagrangian density for the TM1 model including both nucleons and hyperons is written as [3;E3 



B 



+ tK^) 2 - jP^P^ + \ m l^Pn + J2^b (9obo + g^B^i^ + gpp^TB) i> B , (i) 

B 

where the sum runs over all the baryons B—p, n, A, £ 0,± , E 0,_ . In the above Lagrangian a, ui and p are 
respectively the iso-scalar scalar <j, the iso-scalar vector ui and the isovector vector p meson fields. In 
cq. (fTJ) ui^ and p^ v denote, respectively, the field tensors for the uj and p meson fields. 

This model is good enough to describe nucleonic matter and the nuclear saturation point. But it 
is insufficient for the hyperonic matter, because the model does not reproduce the observed strong AA 
attraction. This defect can be remedied by adding two new meson fields with hidden strangeness, namely, 
the iso-scalar scalar a* and the iso- vector vector cj>, which couple to hyperons only [l^]. These fields can 
be identified with the /o(975) and 0(1020) mesons. The corresponding Lagrangian is given by 

Y 

where index Y runs over hyperons only. 

For a complete description of the beta equilibrated cold matter the model should include leptons; 
namely electrons and muons. The leptonic part of the lagrangian is 

C l = ^i(ip-m B )i>i (2) 

The Lagrangian. density of the complete model, which we call TM1YY, is written as 

C TM1YY = C + C YY + C l (3) 

The nucleon coupling constants are chosen from the fit of the finite nuclei properties. The vector 
coupling constants of the hyperons are chosen according to the SU(6) symmetry and the hyperonic scalar 
coupling constants arc chosen to reproduce the measured values of the corresponding optical potentials. 
Below we use the set of model parameters obtained in ref . [17[ • 

We calculate the energy density and pressure for the TM1YY model in the mean field approximation. 
They are given by the following expressions 

H 1 22,^3 I C 4 | 1 2 *2.1 2 2, 3d 4 

e = 2 m ^ <T + 3 CT +2 m ' 7 * <T + 2 TO " W ° + T W ° 
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\™l<& + E B I ' dkk ^ k2 + m B 2 , (4) 



iff 1 2 2 ^3 c 4 1 2 *2 . 1 2 2 , ^ 4 



imXo + + E / ' dk J==== 2 ' ^ 



where m* B = tub — g a BO ' — g a * b&* is the effective mass, vb is the degeneracy factor and fcj? = yj h b — m* B 
is the Fermi momentum of the baryon species B. 



B. Quark Phase 

At higher densities baryons begin to overlap and loose their individuality. In order to describe the 
medium the quark degrees of freedom need to be included. The density inside a compact star is high 
enough to encompass these degrees of freedom. In order to describe the quark phase we adopt the simple 
MIT Bag model [14|, with three flavours (u, d and s). We also add electrons and muons to describe the 
beta equilibrated matter as in the case of the hadronic phase. For the bag model the energy density and 
pressure can be written as 

a 







E ^2 /_ *fcV w / +fc2+B ' ( 6 ) 

f=u,d,s 



pQ - E $ dk-r4=-B, (7) 



where kp = y p 2 f — TWf is the Fermi momentum of quarks with flavor /. For each flavor we choose the 
degeneracy factor vj = 2 (spin) x 3 (color) = 6 and take the following values of quark masses: m u = 5 
MeV, m d = 10 MeV and m s = 150 MeV. 

It is worth noting that because of the negative vacuum pressure (-B in eq.(7)) the Bag model EOS 
always has a zero pressure at a finite baryon density, p* B ~ S 3 / 4 . By this reason the equilibrium 
configurations of the strange quark matter (SQM) may exist even without gravity [18]. They should 
have a sharp boundary with the density jump from p* B to zero. The EOS derived from the NJL model 



has the similar property 19(. 



C. Construction of the mixed phase and EOS 

We are studying the electrically neutral stellar matter in beta equilibrium. Under such conditions the 
chemical potential of a particle species i can be written as 

Hi = Bills + QiVQ (8) 

where Bi is the baryon number of the species i , Qi denotes its charge in units of the electron charge, \xb 
and hq are the baryonic and electric chemical potentials, respectively. Here we assume that neutrinos 
can freely escape from the star. The above equation signifies that only those reactions are allowed which 
conserve charge and baryon number, however strangeness is not conserved. Two independent chemical 
potentials, hb and HQi are found by fixing the baryon and electric charge densities: 

Pb = E B *P l ' PQ = E ® iPl ' ( 9 ) 

i i 
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where pi is the number density of the particle species i . We require electrical neutrality of a star on a 
macroscopic scale, i.e. pq = 0. According to Eq. ([5]). the baryon chemical potential ps equals the neutron 
chemical potential p n and pq is equal to the electron chemical potential p e . At given ps and Pq, the 
quark chemical potentials are found by using the formulae p u = (pb~ 2/iQ)/3 and pd = p s = (pb+Pq) /3 . 

As indicated by the model calculations (see e.g. ref.[20() the deconfinement phase transition, at high 
densities, is of first order in nature. So this transition should produce a MP between a pure hadronic and 
a pure quark phase. There are two ways by which one can construct the MP: the Maxwell construction 
(MC) and the Gibbs construction (GC). Below we consider both possibilities. 

The Gibbs conditions for the mixed phase are 

Pi(vb,Pq) = P 2 (pb-,Pq), (10) 
Pb = Pbi = PB2 , (11) 
MQ = MQi = MQ2 ■ (12) 

Here and below 1 stands for the hadronic phase and 2 stands for the quark phase. 

It is well known that in the case of two chemical potentials the Gibbs conditions (10) - (12) can be 
fulfilled onl y if the coexisting phases have opposite electric charges and the condition of global neutrality 
is imposed [2l|. This condition can be written as 

Pq = (1 - X)p q1 (pb,Pq) + \p Q 2(pb,Pq) = 0. (13) 

Then the volume averaged energy density in the MP is calculated as 

e= (1 - \)e 1 (p B ,p Q ) + \e 2 (p B ,PQ) , (14) 

where A = V2/V is the volume fraction of quark phase. 

Thus, the mixed phase is a very inhomogeneous state of matter with intermittent domains of two 
different phases. Therefore realistic approaches must take into account not only differences in the bulk 
properties of these phases, but also additional contributions to the thermodynamic potential coming from 
the interface energy and electrostatic energy associated with theses domains. First attempt to perform 
such calculations have been done in 12 2j but due to significant uncertainties in the model parameters, as 



< .._</. interface energy, the results are not conclusive yet. As pointed out in ref. it may happen that 
the GC mixed phase is energetically too expensive and may be expelled from the star at all. Then the 
situation is closer to the MC case, where two pure phases are in direct contact with each other. This 
situation corresponds to the Maxwell construction of the mixed phase defined by the conditions : 

Pi(Pb-.Pq) = P2(pb,Pq) (15) 
Pb = Pbi = PB2 (16) 

They mean that the baryon chemical potential is continuous, but the electric chemical potential pq jumps 
at the interface between the two phases. Also, contrary to the Gibbs construction, where the pressure in 
the mixed phase increases with baryon density, Maxwell construction corresponds to constant pressure 
in the density interval of the mixed phase. 

Figures la and lb show the equations of state obtained with GC and MC for two cases, B 1 / 4 = 
180 MeV and 185 MeV, respectively. As can be seen from Fig. la, for MC the mixed phase starts at 
Pb = 0.41 fm -3 and ends at 0.68 fm~ 3 , whereas for the GC these values are 0.23 fm -3 and 0.89 fm -3 
respectively. So in the case of GC the phase transition starts early and the width of the MP region is 
much boader compared to that in MC. As one increases the Bag constant B the width of the MP region 
further increases for GC, as can be seen in Fig. lb. In this case the MP starts at ps — 0.26 fm -3 and 
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FIG. 1: Equations of state for a) B 1/4 = 180 MeV and b) B 1/4 =185 MeV, for Maxwell and Gibbs construction 
of the mixed phase. 



ends at 1.06 fm~ 3 . On the other hand, for the MC the width of the MP region decreases. For this case the 
MP starts at ps — 0.67 fm -3 and ends at 0.83 fm~ 3 . These results will certainly affect the star properties 
which we will discuss in the next section. 

In Fig. 2 we have plotted the particle abundances for all the cases discussed above. In Figs. 2a and 
2b the results are for B 1 / 4 = 180 MeV, for the MC and GC, respectively. From Fig. 2b, i.e. for GC, we 
see that only one hyperon, i.e. A, is present in the medium. On the other hand, for the case of MC (Fig. 
2a) both E~ and A hyperons appear in the medium. This happens because in the case of GC the phase 
transition starts early, and as a result, hyperons can appear only in the MP. However due to the presence 
of the strange quark, hyperon production in the MP is suppressed and as a result £~ does not appear at 
all. For MC the phase transition starts much later allowing the £~ to appear in the hadronic phase. It 
is interesting that the case of B 1 / 4 = 185 MeV exhibits a completely different picture as shown in Figs. 
2c and 2d. Firstly, in the case of GC, the MP region is broader as compared with the case of B 1 / 4 = 180 
MeV. This allows almost all the hyperons, except S°, to be present in the matter, the E~ appears after 
A. For the case of MC, as the MP begins at a higher density, almost all the hyperons are present in the 
hadronic phase. But, contrary to the case of GC, the S~ appears before A. So the particle cocktail is 
rather different for the two constructions. 



III. PROPERTIES OF COMPACT STARS 

Having obtained the EOSs and the particle abundances we now calculate the properties of compact 
stars with the deconfinement phase transition. We treat the matter to be an ideal fluid and obtain the 
star structure by solving the TOV equations with the corresponding EOS as an input (see details in 
[23!]). The masses and radii of stars are calculated as a function of the central baryon density. Results 
of our calculations for MC and GC are shown in Figs. 3a and 3b, for B 1 / 4 = 180 MeV and 185 MeV, 
respectively. 

In Fig. 3a one can see that the maximum mass for the MC is noticeably higher than that for the GC 
case. The values are 1.493M0 and 1.396M0 respectively. This happens because in the MC the phase 
transition starts quite late compared to the GC, that allows the star to stay longer in the hadronic phase. 
Moreover, for the MC case there is a range of central baryon densities, from 0.41 fm -3 to 0.68 fm -3 , 
corresponding to the MP, which are not allowed in the star. As a result there appears a plateau over 
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FIG. 2: Particle compositions for a) B 1/4 = 180 MeV MC, b) B 1/4 = 180 MeV GC, c) B 1/4 = 185 MeV MC and 
d) B 1/4 = 185 MeV GC. 



this range of central baryon densities. As calculations show, stable star configurations are located on the 
left side of the plateau i.e. till the point A. The region from A to B is not accessible inside the star. On 
the right side of the plateau there is a region, from B to C, where the mass decreases with increasing 
central baryon density, that corresponds to be the unstable configurations. Another interesting feature 
of the MC stars is the appearance of the growing portion of the curve between C and D. This gives rise 
to a new family of stable stars (twin stars) around the mass of I.SAMq similar to the situation studied 
in ref. [2o| . However there is no such a stable solution for GC. The main difference of twin stars from 
the normal neutron stars is that they contain a large quark matter core (see also Fig. 5a). 

For B 1 / 4 = 185 MeV the maximum masses are 1.567M and 1.466M for the MC and GC cases 
respectively. The forbidden range of baryon densities for the MC case is much smaller, ranging from 
0.67 fm -3 to 0.83 fm~ 3 only. However, there is no stable twin star solution for this value of the Bag 
pressure. This situation can be explained from the well known fact that a significant quark core can 
appear only if the density jump is large enough (^J. 

The features discussed above are even more obvious if one looks at the mass-radius plots shown in Fig. 
4. As one can see in Figs. 4a and 4b the mass-radius relations are very different for MC and GC stars. 
Furthermore, the appearance of two maxima associated with the twin stars is very obvious in Fig. 4a. 
The cusps in the MC curves correspond to the plateau regions of Figs. 3a, b. 

We have calculated also the baryon density profiles of MC and GC stars for the same central density. 
They are given in Figs. 5a and 5b. Figure 5a shows a sudden jump in the baryon density for the MC 
case whereas in the GC case the profile is smooth. As Fig. 5b shows there is no jump for the case B 1 / 4 
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FIG. 4: Mass-Radius relation of the stars for a) B 1/4 = 180 MeV and b) B 1/4 = 185 MeV. 



= 185 MeV and there is no stable star with a quark core. 



IV. JUMP IN ELECTRON CHEMICAL POTENTIAL 



In this section we discuss possible implications of the jump in the baryon density of the MC stars as 
illustrated in Fig. 5a. Figure 6 shows the radial profiles of the electron chemical potential /i e for the case 
of B 1 / 4 = 180 MeV. As expected, in the GC case fx e evolves smoothly over the mixed phase reaching 
almost zero in the pure quark phase. But in the MC case fi e has a jump at the transition between 
hadronic and quark phases. In the quark core /i e is very low about 10 MeV and is almost constant, but 
in the hadronic phase /j, e jumps to a high value, about 196 MeV. Then at larger radii /j, e evolves slowly 
and finally decrease to small values as one approaches the crust. 

This is an interesting situation in the sense that such a discontinuity in the chemical potential of 
electrons would lead to the flow of electrons across the discontinuity surface from the region with higher 
[i e to the region with lower \x e inside the star. This flow will be terminated by the electric field generated 
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FIG. 5: Comparison of baryon density profiles for MC and GC stars calculated for a) _B 1//4 = 180 MeV and b) 
B 1/4 = 185 MeV. 

due to the charge separation. The equilibrium condition can be expressed as 

Hei - e$i = /i e2 - e$ 2 (17) 

where fiei and \i e i are the electron chemical potentials in the hadronic and quark phases, and $1.2 are 
the electrostatic potentials far away from the discontinuity. 
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FIG. 6: Profiles of electron chemical potential for MC and GC stars for B 1/4 = 180 MeV. 



The profile of the electrostatic potential $(r) over the discontinuity surface can be found from the Pois- 
son equation in combination with the Thomas- Fermi approximation for the electron density p e (r) — ^f^n 
where kp(r) is the local Fermi momentum of the electrons (see details in ref. 25j). In the ultrarelativistic 
limit hi » rn e , /12 >> rn e , the result can be obtained in the analytic form and the maximum electric 
field is expressed as 

^ = (J-elHe2_ x Hel He2 (lg) 
eir Hel + Me2 

For the particular case shown in Fig. 6 we have fi e i = 10 MeV and He2 — 196 MeV, that gives Eq k 3 
MV/fm. This is a very strong field, about 1000 times the critical field needed for the spontaneous 
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electron-positron pair production in vacuum. In fact, this pair production process is Pauli blocked in the 
considered case. Nevertheless, such a strong electric field may lead to interesting phenomena such as e.g. 
generation of strong magnetic field in the rotating stars. The generation of strong electric fields at the 
bare boundary of a quark star was first discussed in ref. 26 1. 



V. SUMMARY 

To summarise, we have compared the Gibbs and Maxwell constructions of the mixed phase in the 
context of deconfinement phase transition in compact stars. For this purpose we have used a RMF model 
(TM1YY) for the hadronic phase and the MIT Bag model for the quark phase. We have found that the 
EOSs are very different in these two cases: the MP region occupies a much broader density interval for 
GC as compared with the MC one. As we increase the Bag constant the width of the MP region in case 
of GC increases whereas it reduces for the MC case. The particle compositions are also found to be very 
different for the two cases. We then use these EOSs to calculate the star characteristics. The maximum 
mass is found to be different for the MC and GC cases. Furthermore, for i? 1 / 4 =180 MeV a stable solution 
corresponding to second family of compact stars is obtained for MC. The baryon density profiles show a 
sharp jump for the MC case. 

We have also studied the behaviour of the electron chemical potential fi e across the star and found 
that it jumps sharply, for MC, at the phase transition boundary. The jump is about 185 MeV for B 1 / 4 = 
180 MeV. We point out that this jump will lead to redistribution of electrons and generation of a strong 
electric field at the phase transition boundary. We are planning to study this interesting effect in the 
future. 
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